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Summary 0- / f '? f 

A theory of crater formation in solids by impact of ultrahigh- 
speed particles is studied from the standpoint of radially sym- 
metric advancing shock fronts. The equations of motion lead to 
a solution based on progressing waves, which leads to a 2/5-power 
law for penetration versus velocity. Some particular results for 
steel agree with Charters’ and Summers’ data. The variations 
of pressure and density with radius are also obtained. 
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pressure 

radius 

radial velocity of impact of projectile 
shock- wave velocity, or “sound” velocity 
particle velocity 

volume, Vo ~ initial volume, V/Vo — relative 
volume 

density (p p , of projectile; p t , of target) 

initial (uncompressed) density 

mass 

adiabatic exponent 
progressing-wave exponents 
progressing-wave parameters 
progressing-wave functions 
V(y - l)/(y + 1 D 
energy 
momentum 
diameter of projectile 


(1) Introduction 


T he danger of ultrahigh- speed impacts of meteoric 
particles on space vehicles has led to an interest in 
theories of penetration of small pellets into metal 
plates. The phenomena; and hence their analysis, are 
complex because different effects predominate in dif- 
ferent parts of the crater. Thus in a zone ahead of the 
projectile, the target material is probably being tri- 
axially compressed. This zone is bounded in front by 
an advancing shock wave and gradually shades off at 
the sides into another zone where the material, in some 
liquid or plastic state, flows sideways along the crater. 
Ultimately, as the angle of deviation from the direction 
of impact is increased, the material has undergone per- 
manent plastic deformation and hills up to form the lip 
of the crater. These zones are shown in Fig. la. The 
transition between the two zones is not meant to be ex- 
act in such a schematic picture. 

Our approach here will be to consider a simplified 
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hemispherical zone, as shown in Fig. lb, and -which, in 
effect, extends the shock front all the way around to 
6 = 90°, with spherical symmetry and radial particle 
motion. 

The direct aim of any penetration theory is, of course, 
the prediction of the form of the crater. In the spheri- 
cal-zone problem this is done by analyzing the state 
of the material ahead of it, leading to a value for the 
diameter and volume, for comparison with experimental 
data. Corrections can then be made for the flow effects 
in Zone II. 

Craters have been produced in copper by high- 
velocity microparticles by Anderson, Doran, Hempy, 
and Kells. 6 Important cratering experiments have been 
performed by Charters and Summers, 3 ' 4 by Partridge, 
Morris, and Fullmer, 5 and by Atkins. 9 Analyses have 
been made by Gehring, 7 and theories of crater forma- 
tion have been reviewed by Allison. 11 Recently, Band 13 
has treated shock propagatio»-in viscoelastic mediums. 

(2) Spherical Shock Wave — Distribution of 
Effects in Time 

Just as the study of the problem is conveniently 
divided into zones, we can divide the sequence of events 
in the crater-formatioti process as follows for detailed 
analysis. 

(1) Initial Stage — here the projectile at impact be- 
comes imbedded just inside the target material and 
generates a region of highly compressed material. The 
problem is to determine the nature of this zone. 

(2) Expansion Stage — the compressed material, act- 
ing in a manner similar to an explosion, expands further 
into the target, generating shock fronts, and forming the 
crater. 

(3) Final Stage— the shock wave decays, permanent 
deformation at the crater stops, and secondary defor- 
mations appear at the back or other parts of the target. 
The main problem is to determine when this phase 
begins. 

Of course, it must be understood that these phases are 
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Fig. la. Region 
around impact cavity. 
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not distinct events, nor may it even be possible in a 
given case to make the separation. 

(3) The Initial Stage 

We lack direct data on the initial stage of crater for- 
mation. Headington and Jaunzemis 17 have analyzed 
this stage upon the basis of conservation laws and the 
known equations of state of the materials. Such a de- 
tailed study provides the initial and boundary condi- 
tions required for Stage 2. However, we shall here make 
simpler assumptions about the initial phase during 
which the projectile becomes imbedded in -the target 
material : (1) the impact is so rapid that there has not 
been time for flow to take place; (2) the time required 
for the pressure to reach its maximum (compression 
time) is negligible, which ’means that Ve can treat the 
problem essentially as one of an explosion at the point of 
impact, 0; (3) no heat or energy is lost during the com- 
pression stage — i.e., the process is adiabatic; (4) the 
pre-impact shape of the projectile is unimportant. 

These assumptions lead us to the so-called ‘‘ballistic” 
model of impact initiation. 3 In Ref. 3 it is supposed 
that immediately on impact the projectile and some of 
the target material are converted to a fluid shell under 
high pressure which propagates radially. In this 
analysis, however, we shall remove the assumption that 
the shell is uniform — part of our effort will be to deter- 
mine the pressure profile and velocity distribution in it 
as a function of radius. 


(4) The Theory of Progressing Shock Waves 


We have previously (Fig. 1) considered the different 
angular zones in cratering. If we restrict the problem 
to spherically symmetric forces, and hence to radial 
velocities, we simplify it to the point at which it can be 
analyzed theoretically. However, since edge effects are 
being neglected, the results will be limited to some 
angle less than 90° to the normal. 

Our model is related to that of a spherical blast wave 
in a gas, and we may make the following assumptions 
about the medium during propagation of the shock 
wave: (1) thermodynamic equilibrium holds (Ref. 2, 
p. 2) — i.e., changes of state are adiabatic, or entropy is 
constant along a particle path; (2) the medium is a 
perfect fluid — i.e., any rigidity effects are neglected; 
(3) the effects of entropy changes are negligible, which 
implies that the pressure is a function of the density 
alone (medium barotropic); (4) there are two alterna- 


tive possibilities: (a) the total energy available for the 
motion is given, or (b) the total momentum available 
for the motion is given. 

(5) The Basic Equations 

The conservative laws for an element of material may 
be expressed in Lagrangian form — i.e., along the particle 
paths — as 
Momentum: 


Mass: 


du/dt — — dp/pdr 


( 1 ) 


dp/dt — —pdu/'dr — 2 pp/r (2) 

equation of state: 

f{p> P) = 0 (3) 

If the medium is assumed to be polytropic with the 
adiabatic exponent 7, we have for Eq. (3) 

f(P, p) = pp' y = const. = A (4) 

These equations, in Eulerian form, become 

u t + uu T -f- (1/ p)Pt ~ 0 (5a) 

Pt + up r + pu T + (2 up/r) — 0 (5b) 

(pp- y )< + u(pp~ y ) T = 0 (5c) 


Eq. (5c) expresses the fact that the entropy is constant 
along the path of an element, although it does not imply 
that the entropy is constant throughout, as in the plane 
case. Instead we have a condition of constant total 
energy. We also note that the only difference from the 
equations for one-dimensional flow is the additional 
term 2 up/r occurring in Eq. (5b), which stands es- 
sentially for the spherical attenuation of the wave. 
This term, of course, is very important in the problem. 


(6) Progressing-Wave Solutions 

We shall use a known general method which reduces 
these equations to a succession of ordinary differential 
equations. By assuming a specific form for the equa- 
tions a set of particular solutions depending on one 
variable will be obtained. These are called ‘‘progress- 
ing waves.” For general details, see Ref. 1, pp. 419- 
433. 

The “progressing wave” solutions are defined to be of 
the form 

u = t e mk) ) (6a) 

P = t s D( |) > with £ = rt~ a (6b) 

p/p = * e em) (6c) 

where a, /3, 5, and e are parameters, and U, D, and P are 
functions to be determined. By introducing the vari- 
able £ we have defined geometrically a family of sur- 
faces £ = const, in the r.bplane, which will play an im- 
portant role in the analysis. Although these are not the 
trajectories of the particles of the medium, we shall see 
that the shock front belongs to this family of surfaces. 
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Fig. 2. Graphical equation of state for iron under high pressure 
(data from Walsh, Rice. McQueen, Yarger). 

We now explore these solutions mathematically by 
substituting Eqs. (6a) to (6c) into the equations of 
motion (5a) to (5c). It becomes possible to eliminate 
the time t wherever it appears by choosing 

^ = 2/3; 0 — a 1 (7) 

We are then led to the set of ordinary differential 
equations 

tD'/D = -(5 + £t/' + 3 U)/(U - a) (8) 

££/' = [-£/([/- a)(U- 1) + 

(6 + 20 + 3Uy)P]/[(U - a) 2 - yP] (9a) 

tP'/P = {(7 - 1 )U(U - 1) - (U - a) 

[(3 7 - l)tf - 2] + [(2/J - yS + d)/ 

{U- a) +2 y]P}/\{U - a)* - yP) (9b) 

where the prime denotes differentiation with respect to 
the single independent variable £. Dividing Eq. (9a) 
by Eq. (9b) gives 

dP/dU = £P'/££/' = F(P, U) (10a) 
where, after simplification, 

F(P, U) = 

P[yU(3ot -1-2 U) +(3 - a)U - 2a + 

(j [2P - (7 - m/(U - a)} + 2 y)P) 4- 

[— U(U - ac)(U - 1) + (8 + 2/3 + 3Uy)P] (10b) 

This is the basic differential equation for progressing 
waves. After the appropriate solution has been found 
for P = P(U), the function £ = £(£/) is found by a 
quadrature of Eq. (9a), and the density function L>(£) 
from Eq. (8). 


These progressing wave solutions, as we shall see, 
provide a sufficiently general mathematical description 
of an expanding cavity which is reasonably consistent 
with the given conditions of initiation of the process. 
There remains the problem of choosing the two parame- 
ters a and 5. 

(7) Boundary Conditions at Shock Front 

We shall narrow down the number of parameters by 
examining the compatibility of our solution with the 
basic Rankine-Hugoniot conditions across a shock front. 
For an undisturbed medium at rest and with U \ , Pi, pi in 
the disturbed medium behind the shock moving with 
velocity c, these relations are (Ref. 1, pp. 123-4). 

Pi(c — ui) — pqC = 0 (11a) 

p,«i(c — «i) — pi = 0 (lib) 

Pi ( \ uF H ( c - Mi) - piU ! = 0 (11c) 

\2 7—1 pij 

Since £ = rt~ a = £i = constant along the shock front, 
c — dr/dt — art~ l — From Eqs. (6a) to (6c) 

and (7), the shock conditions (11a) to (11c) become 

(A) t* + %D(a - U) - p 0 a£/ = 0 

(B) < H2 ^ 2 Dl r (a - U) - t 8+2 %*DP = 0 

(C) + P/(y ~ 1)} X 

&(a - U) - t 5 +*%*DPU = 0 

We note that the time factor cancels out of (B) and (C), 
but to secure independence of time in (A), it is neces- 
sary to make 

8 = 0 (12) 

With this condition, and Eq. (8), the assumed form for 
the progressing-wave solutions reduces to 

h = <r/f)£7(& p = (r/t)*D(QP(Q) 
p = L(£), p/p = (r//) 2 P(£) > (13) 

with £ = rt " ) 

This solution shows that on the shock front or free 
surface, where £ is constant, the physical quantities 
such as velocity, pressure, density, and wave velocity 
are constant on the rays r/i — constant. Also, this 
correctly makes the functions (13) dimensional, for, 
with £ having the dimensions of length per (time)", and 

£ = [LT~ a ] D(Q = [ML- 3 ] density 
U = [1] P(£) = [1] 

u y py p, p/p correctly are found to have the dimensions 
of velocity, pressure (. ML~ l T~ 2 ), density, and velocity 
squared, respectively. 

(8) Energy Condition 

The assumption of constant total energy for the 
cavity expansion process is a reasonable one for a fast 
process/' prbviij&i } certain secondary effects are neg- 
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lected. With £ = £1 representing the shock front at a 
time t } the total energy in the fluid shell (kinetic + po- 
tential) at time t is given by 

E(t) = | { l / 2 PU 2 + [p/(y — 1) ]} 2tt r 2 dr (14) 

Jro = frl a 

where r 0 is the inner radius of the shell (Fig. lb). The 
second term represents the work done in compressing 
the material, based on the assumed poly tropic relation 
p = Ap y . Using the substitutions in Eqs. ( 6 a) to ( 6 c), 
with t constant, we obtain 

E(t) = 2 TTt s+5a ~ 2 f* Y- £7* + — — p)i*Dd( (15) 
[2 7“1 / 

Since the integral is independent of t , we make the 
energy independent of time by satisfying the relation 

5 + 5a - 2 = 0 or 5 = 2- 5a (16) 

The complete set of exponents is now 

ol = 2 U <= = ~ 6 A) Q7) 

P = - 3 A 5=0 j {U) 

(9) Momentum Condition 

If it be supposed that the momentum of the expand- 
ing cavity material remains constant, we obtain a 
different set of values for the exponents. The axial 
component of momentum of the moving material is 
given by 

M = f V P ucos6dV= ? D Udt 

J V o J £o 

so that, with 5 = 0, time-independence requires that 

4a - 1 = 0 


or 




The question of energy versus momentum is a very im- 
portant one in the current literature. Important ex- 
periments have recently been conducted by Culp 8 and 
Feldman , 10 and the problem has also been discussed by 

Lavrent’yev . 12 


(10) Solution for Impact on Steel 

We shall work out the results for the progressing- 
wave theory for a typical material — iron — for which 
some convenient data are available. 

Equation of State for Iron 

Fig. 2 shows the equation of state for iron . 2 A log-log 
plot, Fig. 3, shows two ranges of pressures in which a 
polytropic law of the form of Eq. (4) can be fitted. 
These are, for the high- and intermediate -pressure 
ranges, respectively, 

pp- 9 = 6.9 X 10 - 7 , 0.85 < V/Vo < 0.93) ( . 

pp~ 16 = 1.3 X 10 - 13 , 0.76 < V/Vo < 0.85 j { U) 

where p is measured in kilo bars and p in gm/crn 3 . Be- 
low about 80 kilobars we have a transition to elastic or 
elastic-plastic behavior. The nature of this transition 
is considerably uncertain and will not be discussed 
further here. 

Unlike gases, the high value of y — i.e., the compara- 
tively small changes in density — shows up in that the 
strong shocks which occur in the case of gases, when 
the density changes by a factor of about 6 , are not 
present. 


(11) P U Diagram — Intermediate Pressure 
Range — Constant Energy 

If we assume constant energy in the cavity-formation 
process, the 2 / 5 -power law holds, and with the numeri- 
cal values given by Eq. (17) the differential equation 
for P [Eq. (10a)] becomes, for 7 = 16, 

dP/dU = [N(U) + Q(U)P]/[R(U) + S(U)P]} 

where 

N(U) = —32 U 2 + 5.8Z7 - 0 . 8 ; 

R(U) = -U(U - 2/5) (Z7 - 1) 

Q(U) = 32 - [6/(517 - 2)]; S(U) = (4 SU - 6/5) 

A family of integral curves is sketched in Fig. 4, made 
from slopes calculated on a grid network interval of 
0.05. Certain portions have been plotted at finer in- 
tervals for a more accurate determination of the initial 
behavior of this solution curve. 

One type of solution is obtained by starting from 
the shock front and applying the known shock transi- 
tions for pressure, density, and velocity. 

If we consider a small surface element of the shock 
front, we can neglect the sphericity in its immediate 
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Fig. 4. Integral curves on PU-diagram for progressing waves. 


neighborhood and obtain the transition relations giving 
pressure, density, and velocities immediately behind it. 
From the Rankine-Hugoniot equations we obtain 


Pi/po = (7 + 1)/(y — 1 ) = i/m 2 
pi = PoC 2 (l — M 2 ) 

Ml = c( 1 - M 2 ) \ 

[dp /dp]' 1 ' = m(1 + M 2 )' /! cj 


( 20 ) 

( 21 ) 


The last quantity is conveniently referred to as the 
“sound speed” in shock-wave studies. 

Since £ = £1 on the shock front, we have, just behind 
it [referring to Eq. (17) and c = dr/dt = <x(r/t)] 


U(id = «( 1 - M 2 ) ( 22 a) 

^(£0 = po/m 2 (22b) 

P(£i) = - P 2 ) ( 22 c) 


with a = 2/5 or 1/4, according to whether the energy or 
momentum condition holds. 

The right-hand sides of Eqs. (22a) to (22c) give us, 
for a specified material, a definite point in the P-U 
plane, through which a single solution curve is deter- 
mined in general. We may refer to this point as our 
“initial point,” and proceed to draw the solution curve. 
Note that the constant £ 1 , still undetermined, is not 
needed for this. 

For iron under pressure with 7 = 16, /x = 0.939, and 
C/(£ 1 ) = 0.04758, P(£i) = 0.016609 


(12) Determination of £ and D Functions 

Having determined the relation P(U) as given by the 


solution curve on the P-U diagram, we have effected 
one quadrature. We now perform a second quadrature 
by referring to the differential Eqs. ( 8 ) and (9a). In 
both of these the variables are separated, so a straight- 
forward procedure leads directly to the £ and D func- 
tions, shown graphically in Fig. 5. 

Note that £ is determined to within the multiplicative 
constant £ 1 , its value at the shock front. From the dif- 
ferential equation for the density p = D(£), we see that 
p also admits of a free multiplicative constant p = pi 
just behind the shock front. This value is, of course, 
specified by Eq. (11b), which, for iron (po = 7.84 
gm/cm 3 ) is 

Pi = 8.9 gm/cm 3 

We see that the density of the material decreases after 
the shock front has passed, at first rapidly, then it 
levels off until U = 0.341, where it falls so rapidly to 
zero as to create almost a discontinuity there. The 
same occurs with £. 

The interpretation of the curves and their ultimate 
jump is very important to our problem. Since £ de- 
creases with radius r when time is held constant, the 
curves show that the density drops behind the shock 
front, very rapidly at first, while the velocity increases. 
Ultimately there is a sudden density drop to zero, indi- 
cating a rather sharp boundary from the material to 
free space. In Fig. 5 this is at U — 0.341, which would 
represent the crater boundary on the basis of the theory 
worked out here from the polytropic law. Unfor- 
tunately we can not use this discontinuity, as the com- 
puted density of the material has by this time dropped 
considerably from its free-space value po. Since there 
is no evidence that any solid in these processes ever 
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Fig, 5. Profiles of density and £ parameters. 


has a density less than p 0 , we must find a way to de- 
viate from the polytropic law at low pressures. The 
simplest way is to level off the density at p = p 0 . The 
effect of this will be examined in the next section. 

(13) The Particle Trajectories 

In order to study the motion of the crater surface, 
there remains only one more construction to be made, 
that of the r, 2-plane, showing the particle trajectories. 
(By “particle” we mean actually all the particles on a 
surface of radius r .) This is done by linear elements, 
computed from the relation 

dr/dt = u(r , t) = (r/2)t/(£), £ = r2~ 2/5 (23) 

which is Eq. (13). The result is shown, for our numeri- 
cal data, in Fig. 6. It is convenient, because of the 
rjt factor, to plot the linear elements along rays r/t = 
const., as shown. 

Of particular interest in this diagram is the shock- 
front path, given by 

r = 

with £i yet to be determined. This is done from the 
energy condition (see below). 

(14) The Pressure 

The pressure variation in the medium is given by [see 
Eq. (13)] 

P = <r/t)*D{QP(d (24) 

which shows the same discontinuity at the crater sur- 
face as the density function does. The values [Eq. 
(24)] just inside the medium give a nonzero pressure 
variation 

P = O'/ 2) 2 (0.0004) Pl £i 2 

as shown by the dotted line in Fig. 7. This violates the 
constant-energy condition (since some work is con- 
tinually being done at the boundary), but it shows ap- 
proximately the impulsive-type pressure which the 
conditions of the problem require. On the other hand, 



Fig. 6. Motion of fluid shell and crater surface on rt-diagram 
for impact at 20,000 fps. 


just beyond the jump, the vanishing of D gives us zero 
pressure as the energy condition requires. We may 
properly regard this as a “free” surface. 

The only way to come out exactly with a free surface 
without having the discontinuity in pressure just de- 
scribed, is for the solution curve in Fig. 4, to end at 
U = 0.4, P = 0, as shown by the dotted line in the 
figure. Only then would we be able to say that the 
progressing-wave solutions fit the physical conditions 
required by the problem. As the plotting here is very 
sensitive to the initial conditions, we cannot say now 
for sure whether failure to reach this point is due to in- 
accuracies, or whether the initial conditions need to be 
changed slightly. However, our solution approaches 
the required point closely enough to delineate the crater 
surface behavior, which is hardly sensitive at all to this 
discrepancy. 

Physically, the pressure drops very rapidly with de- 
creasing radius from the shock region. Just how or 
where it drops to zero is uncertain, because we are 
neglecting the departure of the material from the poly- 
tropic law at low pressures. This does not appear to 
matter too much in locating the locus of the crater 
surface. 
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(15) The Energy 

The above analysis of shock-wave propagation has 
been carried out on the basis of constant energy in the 
process. Any estimate of ultimate crater size must de- 
pend on the value of this energy, which unfortunately 
is not and cannot be provided by any theory of propaga- 
tion alone, as it is one of the input parameters in the 
problem. The input energy consists of the impact 
kinetic energy of the projectile. This becomes par- 
titioned among various physical effects during the 
cratering process — impact kinetic energy of projectile: 
(1) binding energy of materials; (2) permanent defor- 
mation of plastic flow; (3) heat; (4) light; (5) splashing- 
out of pieces of projectiles; (6) shock wave. We are, 
of course, interested in the fraction of energy available 
for item (6), but there is little information available on 
this point. Some data are given by Partridge. 5 The 
effect of the energy fraction on the results will be dis- 
cussed in the next section. 

We note that in the thin fluid shell just behind the 
shock front the calculations shown previously furnish 
the following values 

U = 0.047, V 2 U 2 = 0.00111, P/(y - 1) = 0.00111 
so that the kinetic energy of the fluid shell is half that 
of the projectile, in conformity with the ballistic model 
for impact (Ref. 3, p. 11). 


(16) A Representative Cratering Condition 

Consider a steel pellet 3 /ie-in. (0.738 cm) in diameter 
striking a thick steel plate at 20,000 ft per sec. This 
value, which happens to be especially critical for ter- 
minal ballistic effects 14 represents a ratio u/c — 1.025, 
which just exceeds a Mach number of 1. This repre- 
sents, with 

p p — 7.84 gm/cm 3 (iron) 

an impact energy (all kinetic) of 

E = y 2 (0.443 gm)(6.10 X 10 5 cm/sec) 2 - 

8.25 X 10 10 erg 

If we suppose all this energy enters the shock wave (no 
loss basis) the integration of the energy expression 
[Eq. (14) ] for our numerical values gives 

E = 2r(3.80 X 10- 4 ) Pl £i 5 - 8.25 X 10 10 ergs 

With pi = 8.90 gm/cm -3 , we obtain 

= 329 cm sec -2/5 

The pressure function behind the shock front is 

p = t^DP X pi£i 2 = 1.60 X 10' 5 ^ /5 kilobars 

with r in cm and t in sec. Further, for 

t = 1 microsec, p = 254 kb 
t — V 2 microsec, p — 582 kb 

This shows that the high-pressure regime (during which 
the crater is being formed) is over in about 1 microsec. 

Fig. 6 is a plot of radius vs time which shows the 
motion of the shock front and trailing surface, and some 
of the particle paths. We note that after the density has 
fallen to its free-space values, at the trailing front, the 
particles still have a residual forward motion, which is 
effectively over by at most 10 microsec. The crater 
surface is stopped at the value 

r ^ 0.37 X 10-^i cm 

On the basis of no energy loss (£i = 329) we obtain, 
as a maximum crater radius, 

r = 1.2 cm 

This must actually be reduced proportionately to the 
fraction of energy available. 

Following the notation of Charters and Partridge, 

p/d = 1.63 

where p = penetration distance and d = diameter of 
projectile. Referring to Charters’ curves of p/d vs 
V/Cy where V = velocity of projectile and c = sound 
velocity, this point lies on the fluid impact curve, in 
the transition region (Fig. 8). 

The experimental points for the impact of iron into 
iron, as obtained by the University of Utah, do not 
conform to the 2/3 law. However, these points only go 
up to V/c — 0.5. The analysis herein, which provides 
a 2/5 power law, is flatter than the 2/3 law. Which 
law really holds for V/c > 1 we do not yet know. 
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Fig. 8. Comparison of 2/5-power law with experimental points 
from Univ. of Utah. 

(17) Conclusion 

Generally speaking, the progressive -wave analysis 
appears to give promise of reasonable agreement with 
experimental values in spite of the many assumptions 
needed and uncertain values of the parameters re- 
quired. Our aim will be later to extend this method to 
other materials and to carry out an analysis based on 
momentum also. 

For the iron-iron collision at 20,000 ft per sec, the 
crater is formed in about 1 microsec. 

The pressure is a maximum just behind the shock 
front. There is then a trailing wave in which the pres- 
sure drops rapidly. Just at the crater, the material 
has been decompressed. From Fig. 6 we see that the 
shock front starts out at about the impact velocity and 
slows down to the elastic -wave velocity during the first 
microsecond, while the crater surface penetrates at a 
much slower speed. 

An interesting side feature is that if we had allowed 
the pressure release front to stop suddenly at 1 microsec, 
as marked by point A in the figure, we would come out 
with about the same value for the crater radius. 

The release front shown in the figure is somewhat 


idealized, of course, since the transition from the high- 
pressure region to the elastic one is actually gradual. Of 
the many questions concerning the cratering process, 
our immediate aim has been directly toward increasing 
our understanding of how shock waves are propagated 
into a material ahead of a developing crater. 
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